Surface critical behavior of fluids: Lennard-Jones fluid near weakly attractive 

substrate 
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The phase behavior of fluids near weakly attractive substrates is studied by computer simulations 
of the coexistence curve of a Lennard-Jones (LJ) fluid confined in a slitlike pore. The temperature 
dependence of the density profiles of the LJ fiuid is found to be very similar to the behavior of 
water near hydrophobic surfaces (Brovchenko et al. J.Phys.: Cond.Matt. 16, 2004). A universal 
critical behavior of the local order parameter, defined as the difference between the local densities 
of the coexisting liquid and vapor phases at some distance z from the pore walls, Ap(z) = (pi(z) - 
p„(z))/2, is observed in a wide temperature range and found to be consistent with the surface critical 
behavior of the Ising model. Near the surface the dependence of the order parameter on the reduced 
temperature t = (Tc - T)/Tc obeys a scaling law ~ t^'^ with a critical exponent /3i of about 0.8, 
corresponding to the ordinary surface transition. A crossover from bulk-like to surface-like critical 
behavior with increasing temperature occurs, when the correlation length is about half the distance 
to the surface. Relations between the ordinary and normal transitions in Ising systems and the 
surface critical behavior of fiuids are discussed. 

PACS numbers: 05.70.Jk, 05.70.Np, 64.60.Fr, 64.70.Dv 



I. INTRODUCTION 



The presence of a solid boundary effects the properties 
of a fluid, which becomes spatially heterogeneous normal 
to the boundary. Besides density oscillations in the close 
vicinity of the surface due to packing effects, the surface 
produces perturbations which intrude into the bulk fluid 
on a distance scale proportional to the bulk correlation 
length and, therefore, gain increasing importance near 
the critical point even far away from the surface. The 
critical behavior of the local properties of the fluid near 
the surface should follow universal scaling laws, which 
are different, however, from the scaling laws of the bulk 
fluid. On strict theoretical grounds simple scaling equa- 
tions are valid only in the immediate neighborhood of the 
critical point. However, both experimental and computer 
simulation studies of the coexistence curves of bulk fluids 
evidence that contributions from the asymptotic power 
laws remain dominant over a surprisingly wide tempera- 
ture range. The same situation could also be expected for 
the surface critical behavior. In this case the knowledge 
of the surface critical behavior of fluids opens the possibil- 
ity to predict the density profiles and related properties 
of various fluids near solid boundaries in a wide range of 
thermodynamic conditions. 

Since bulk fluids belong to the universality class of the 
3D Ising model, it is natural to map the surface criti- 
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cal behavior of fluids onto the surface universality classes 
of the Ising model. In the presence of a non-zero sur- 
face fleld (hi ^ 0, the case relevant for fluids) in the 
Ising model a wetting transition ultimately occurs at 
some temperature below the bulk critical temperature 
Tc Q IM IE Above the temperature of the wetting 
transition the magnetization of the single phase, remain- 
ing near the surface, is predicted to follow the law of 
the so-called normal transition 's', 'g'I with dominant reg- 
ular behavior, i.e. with a nonzero constant followed by a 
leading linear contribution, proportional to the reduced 
temperature r = (Tc - T)/Tc, which is a measure of the 
temperature deviation from the critical temperature Tc. 
The available computer simulations of Ising systems at 
J^i 7^ Qi S 01 did not allow to analyze quantitatively 
the temperature dependence of the surface magnetization 
both below and above the wetting temperature. 

There are no experimental studies of the temperature 
evolution of the density profile near a solid surface along 
a liquid-vapor coexistence curve, which would allow to 
explore the surface critical behavior of one-component 
fluids. Experimental studies of the order parameter in 
binary liquid mixtures near a solid surface (difference be- 
tween the concentrations of the two coexisting phases) 
evidence that its temperature dependence near the wall 
follows a power law ~ t^'^ with a surface critical expo- 
nent « 0.8 mm. This value is close to the critical 

exponent of the ordinary transition of the Ising model at 
hi = 0, that describes the temperature evolution of the 
magnetization in a surface layer p5 fTlj . 

It is usually supposed that the order parameter of one- 
component fluids and binary fiuid mixtures near a surface 
obeys the laws of the ordinary transition below the wet- 
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ting (drying) temperature 0, 0, 0|. In the case of 
a weakly attractive surface the long-range fluid-wall in- 
teraction suppresses a drying transition ITR , and 
so the order parameter near the surface should follow a 
temperature dependence 
temperature. 



up to the bulk critical 



In fact, such a behavior of the order parameter was 
observed in computer simulation studies of the liquid- 
vapor phase transition of water near hydrophobic sur- 
faces 13 • The temperature evolution of the densities 
Pi^v{z,t) in coexisting liquid (subscript 1) and vapor (sub- 
script v) phases near a substrate could be described by 
the following equation: 



p,,„(z, r) = {Ao + AiT + A2T^...) ± Bit/^ 



(1) 



where Aq = Pc{z) = p{z, r = 0) is the local density at 
the critical point, Ai, A2... and Bi are local system de- 
pendent amplitudes. The asymmetric contribution in the 
bracket (diameter of the coexistence curve) is the same 
in both coexisting phases and may include also a singular 
contribution ^ r 0,!^; with critical exponent a = 

0.109 19,]. Additionally, a singular contribution due to 
pressure mixing could also be expected, as in a bulk fluid, 



where it is 
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l2C|. Note, however, that contrary to 



the bulk case, the exponents of both of these singular 
terms exceed 1 and, therefore, the regular linear term re- 
mains the most important asymptotically. The last term 
in Eq.(l) represents a symmetric contribution to the den- 
sities of the coexisting phases (has opposite signs in two 
phases), and thus describes the order parameter Ap = 
{pi - Pv)!'^! which becomes zero at the critical point at 
any distance from the surface. The value of the critical 
exponent fi\ was found close to the value « 0.8 of the 
ordinary transition |l7l Il8j| . 

The densities of the coexisting liquid and vapor phases 
of water in a layer of molecular width near a hydrophobic 
surface follow Eq.(l) in a wide temperature range: 0.5 < 
T < 0.05 Moreover, the intrusion of the surface crit- 
ical behavior described by Eq.(l) deeper into the bulk 
fluid was found to be governed by the bulk correlation 
length. So far, this was the only computer simulation 
study of the surface critical behavior of the order param- 
eter of a fluid. Naturally, the universality of the surface 
critical behavior, reported for water near a hydrophobic 
surface should be tested for other fluids and, also, 
for various strengths of fluid-wall interaction, below as 
well as above the wetting (drying) temperature. 

In this paper we present a study of the surface critical 
behavior of a Lennard-Jones (LJ) fluids near a weakly at- 
tractive wall by computer simulations of its liquid- vapor 
coexistence curves in a slitlike geometry. The surface 
critical behavior is explored by considering the temper- 
ature evolution of the density profiles of the coexisting 
liquid and vapor phases. 



II. METHOD 

Gibbs ensemble Monte Carlo (GEMC) simulations ^ 
were used to simulate liquid-vapor coexistence curves of 
bulk and confined LJ fluid having interparticle interac- 
tions of the form: 



UL.j{r) = 4e (a/r)'' - {a/rf 



(2) 



where e measures the well depth of the potential, while a 
sets the length scale. The potential was spherically trun- 
cated at a radius 2.5 u and left unshifted. No long-range 
corrections were applied to account for effects of the trun- 
cation. The density p used in the present paper is the 
number density scaled by , while T is the temperature 
scaled by e/kB, where ks is Boltzmann's constant. 

Simulations of the bulk liquid-vapor coexistence were 
performed at 51 temperatures from T = 0.60 to T = 1.17. 
The two lowest temperatures (T = 0.60 and 0.65) were 
below the bulk triple-point temperature of the LJ fluid 
(the values 0.687 0.689 [H and 0.692 Q were re- 
ported in the literature). The total number of molecules 
in liquid and vapor phases N; -I- N^ was about 2000 at 
low and 1500 at high temperatures. The edge size L of 
the cubic simulation box for the bulk liquid was about 
12 a and the number of molecules varied from ~ 1900 
to ~ 750 with increasing temperature. The number of 
molecules in the vapor phase N„ varied from ~ 100 at 
the lowest studied temperature to ~ 750 near the critical 
temperature. As a strong difference between the number 
of particles in the two simulated phases may distort the 
shape of the coexistence curve |25ii, ,26| , at high tempera- 
tures we always tried to keep Nj, close to N;. 

To estimate the effect of the system size on the coex- 
isting densities, at some temperatures the liquid- vapor 
coexistence was simulated for larger systems, where the 
total number of particles in the liquid and vapor phases 
N; -I- N^ was about 7500, which results in a box size L 
for the liquid of about 20 a. To test the effect of differing 
numbers of molecules in the two simulation boxes of the 
coexisting phases j2^ |2^ we repeated the GEMC simu- 
lations at T = 1.14 and T = 1.15 with N, « 3N^, and N; 
« 0.7 Ny, respectively. We found this effect neghgible, 
obviously due to the large system size. 

To study the surface critical behavior, the LJ fluid was 
confined in a slitlike pore with structureless walls. Each 
wall interacts with particles of the fluid via a long-range 
potential comprising a single plane of LJ molecules: 



U^{z) = 4e f 0.4 (cr/z)^° - (a/z) 



(3) 



where z measures the distance to the wall and the pa- 
rameter f determines the strength of the fluid-wall inter- 
action relatively to fluid-fluid interaction. No truncation 
was applied to Uuj(z). In the present paper we report the 
results obtained for a pore of width H = 12 a and / = 
0.3, which corresponds to a weakly attractive (strongly 
solvophobic) surface. Note, that this system is very sim- 
ilar to the one, recently studied in equilibrium with the 
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TABLE I: Densities of the coexisting phases of the bulk LJ 
fluid (interaction potential truncated at 2.5 a, unshifted), ob- 
tained for the system with i ~ 12 cr. 



FIG. 1: Liquid- vapor coexistence curve and its diameter, ob- 
tained for the bulk LJ fluid (interaction potential truncated 
at 2.5 (T and unshifted). Closed circles: box size i ~ 12 cr, 
open circles: L « 20 cr. Open squares, triangles and diamonds 
indicate the data points, reported in Refs.l27j|. |2^ and |29||. 



critical bulk L J fluid |30j . The total number of molecules 
in the liquid and vapor phases N; -I- N„ in the pore varied 
from 2700 molecules at low temperature to 2000 at high 
temperatures. The lateral size L of the simulation cell for 
the liquid was about 17 a and so the ratio L/ff exceeded 
1. To explore the effect of the system size, several tem- 
perature points were simulated also for a larger system 
with N; -I- Ni, of about 8000, that corresponded to L 
34 CT. 

The efficiency of the molecular transfers between the 
simulation cells was improved by early rejection of in- 
sertion attempts which would lead to strong repulsion 
|3lj| . For each temperature point the number of success- 
ful transfers per particle between the coexisting phases 
of the bulk as well as of the confined fluid was about 100 
to 200, while the probability of successful transfers varied 
from 0.5 to 10%. Close to the critical temperature GEMC 
simulations are limited by identity exchanges between the 
simulation boxes. This strongly increases the error bars 
of the estimated densities of the coexisting phases and 
determines the high-temperature limit of the simulated 
liquid-vapor coexistence curve. We do not give the den- 
sities of the coexisting phases, when identity exchanges 
occured (except for the highest-temperature point T = 
1.17 in the bulk LJ fluid with L « 12 a). The density 
profiles of the coexisting liquid and vapor phases and of 
some supercritical states in the pore were obtained by 
subsequent Monte Carlo simulations in the NVT ensem- 
ble. 
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III. RESULTS 



Bulk coexistence curve 



The coexistence curve of the bulk LJ fluid was studied 
by various simulation methods and its critical parameters 
were found to be strongly sensitive to the details of the 
interaction potential (cut-off, long-range corrections, use 
of shiftedpotential) I21I MMjlilHilslllllllIll 
111 S IM El El EI H il- For example, the 
values of the critical temperature Tc, reported for a LJ 
potential, which was truncated at 2.5 a but not shifted 
(1.1876(3) ,28,], 1.186(2) ^4^, 1.186 ^ 1.1879(4) ^Ij), 
are essentially lower than the values, reported for a full 
LJ potential (1.310 [H, 1.3120(7)[4|, 1.3145(2) E3), 
1.3207(4) E3)- The critical density pc is less sensitive 
to the truncation of the pair potential: values of Pr. from 
0.314 to 0.316 were reported for the fuU potential [MEl 
lil . E3 1 whereas for a L J potential , which was truncated 
at 2.5 a and not shifted it varies from 0.316 to 0.320 

IMIMIillii. 
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The densities of the coexisting phases of the bulk LJ 
fluid with truncated (at 2.5 a) and unshifted intermolec- 
ular potential, obtained in the present studies for the 
system with L w 12 cr and L w 20 cr, are plotted in Fig.l. 
In Table I the values are given for systems with L « 12 
a. Our results are in good agreement with the available 
data for the same interaction |28L l29j | , which are also 
shown in Fig.l. A more sensitive comparison is based on 
the dependence of the order parameter Ap = {pi - Pt,)/2 
on the reduced temperature r. This needs knowledge of 
the critical temperature T^,. The value of T^, could be es- 
timated from fits of the order parameter to the extended 
scaling equation 



Ap = BoT'3(l + BaT^ + S2AT 



2A 



(4) 



where (3 = 0.326 is the critical exponent of the bulk co- 
existence curve, A = 0.52 is a correction-to-scaling 
exponent, which accounts to the deviations from the 
asymptotic critical behavior, and are the system- 
dependent amplitudes. Fits of Eq.(4) with the leading 
term only (fits 1 and 2, Table II) provide a rather satisfac- 
tory description of the data. In particular, if the temper- 
ature is fixed at the most accurately estimated value Tj, = 
1.1876 [111, the obtained value of /3 = 0.324 is close to the 
value of the 3D exponent of the Ising system (3 = 0.326. 
Two correction terms should be included into Eq.(4) in 
order to achieve the best fits of our simulations data (fits 
3 an 4, Table II). "Best" fit means, that addition of fur- 
ther correction terms does not improve its accuracy. 




TABLE II: Values of the parameters obtained from fits of 
Eq.(4) to the order parameter Ap. represents the mean- 
square deviations normalized to a dispersion of the simulated 
data of 0.001 with a confidence fimit of 0.95. The fixed pa- 
rameter are shown italic. 





/3 


Tc Bo Ba 


B2A 


2 

X 






bulk LJ fiuid, L « 12 ct 






1 0.324(1) 


1.1876 0.5450(10) 




2.55 


2 


0.326 


1.1873(4) 0.5460(5) 




2.60 


3 


0.326 


1.1881(6) 0.5247(50) 0.216(42) 


-0.247(40) 


0.74 


4 


0.326 


1.1876 0.5298(21) 0.179(20) 


-0.213(22) 


0.74 






bulk LJ fiuid, L « 20 cr 






5 


(0.326) 0.1892(20) 0.524(19) 0.209(70) 


-0.233(82) 


4.55 






confined LJ fluid, H = 12 cr, i ~ 


17 a 




6 


(0.326) 


1.1475(16) 0.322(13) 1.30 (21) 


-0.895(181) 2.04 


7 


(0.125) 


1.1443(25) 0.097(9) 6.8(1.1) 


-3.05(65) 


2.03 



a) are in good agreement with the value Tc = 1.1876(3) 
obtained using a histogram reweightingmethod with sub- 
sequent mixed-field finite-size scaling 5285 . Use of the lat- 
ter value of Tc does not worsen the quality of the fit (see 
fit 4 in Table II). Taking into account, that the method 
used in 28] provides the most reliable estimates of the 
critical parameters for the considered LJ fluid, we accept 
the value of the critical temperature Tc — 1.1876 for the 
subsequent analysis. Note, that a recent improvement of 
the finite-size scaling method allows a more precise de- 
termination of the critical density but does not influence 
noticeably the critical temperature [4^ . 
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FIG. 2: Log-log plot of the order parameter Ap vs reduced 
temperature r for the bulk LJ fiuid with L ~ 12 a (solid 
circles) and L « 20 cr (open circles, shifted vertically). The 
data from [2^ . |2^ and |29|| are also shifted vertically and 
shown by open squares, triangles and diamonds, respectively. 
The solid lines show the asymptotic critical behavior of the 
Ising model with a slope /3 = 0.326 and parameters of fit 2 
(Table II). 

The values of the critical temperature Tc obtained from 
the flts of our data (Tc = 1.1881(6) for the system with L 
« 12 cr and Tc = 1.1892(20) for the system with L « 20 



FIG. 3: The effective critical exponent /3^^^(r) defined by 
Eq.(5) vs reduced temperature r for several model systems: 
Ising model [3 . full LJ p otential |i33i1 with critical tempera- 
ture Tc = 1.312 from square-well fluid ^ and LJ fluid 
with the potential truncated at 2.5 a (this work). 

The temperature dependence of the order parameter 
Ap is shown in Fig. 2 in double logarithmic scale. All 
data points from Ref.j23| and Ref.|23] and most of the 
data points from Ref.'2^ perfectly agree with our data. 
Two high-temperature points, reported in p^ . slightly 
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deviate upwards from the asymptotic critical behavior, 
probably due to the relatively small system size, used in 
■ Note also, that a slight deviation downward of the 
highest temperature point in our set for the system with 
L « 12 cr may be due to the identity exchanges between 
the simulation boxes, which effectively decrease the dif- 
ference between the estimated densities of the coexisting 
phases. 

The order parameter behavior closely follows the 
asymptotic critical behavior with the Ising critical expo- 
nent P — 0.326 in a surprisingly wide temperature range 
(solid lines in Fig. 2). A more detailed analysis could be 
achieved by considering the temperature dependence of 
the effective exponent P^-^-^ , that describes the behavior 
of the order parameter Apir) at a given temperature (or 
temperature interval) and may be defined as 

P''ff{T)^dln{Ap{T))/dlnT. (5) 

P''^^ can be derived by numerical differentiation of the 
Ap(t) or by analytic differentiation of Eq.(4) with the 
appropriate values of the fit parameters. We applied 
the latter method to plot the temperature dependence of 
peff fj^jj^ studied in the present paper, and 

for some other model systems (see Fig. 3). The deviations 
of the truncated LJ fluid from the asymptotic critical be- 
havior (horizontal line) are the smallest among the avail- 
able simulated systems. The full LJ fluid with long-range 
intermolecular potential shows strong positive deviations 
towards the mean- field critical behavior, which is char- 
acterized by the critical exponent /? = 0.50. The Ising 
model and square-well fluid with essentially short-range 
interactions show a strong crossover to a regular behav- 
ior, which corresponds to P"^^^ 0. So, the coexistence 
curve of the truncated LJ fluid shows really a very weak 
crossover both to the mean-field criticality and to the 
regular behavior and, therefore, it may be a suitable and 
promising model to study the surface critical behavior in 
a wide temperature range. 

The diameter of the coexistence curve contains both 
regular contributions (critical density pc, linear term, 
quadratic term, etc.) and singular terms (~ t^~", ~ 
T^^, etc). Taking into account, that the singular term 
~ T'^^ could be very small in simple fluids oO], not ex- 
tremely close to the critical point, the equation for the 
diameter can be written as 

Pd = {pi + P.)/2 = Pc + Bit + Si^^r^-" + .... (6) 

The diameter of the simulated coexistence curve of the 
bulk LJ fluid is shown in Fig. 4 and the various fits of 
the diameter values to Eq.(6) are collected in Table IIL 
The diameter is essentially linear in a wide temperature 
range (see dashed line in Fig. 4 and linear fits 1 and 3 
in Table III) and only a few (4 or 7) points at the high- 
est studied temperatures deviate upwards from the lin- 
ear dependence. These points show clearly scattering 
around the average values, which give pc ~ 0.335. The 
fits of the Eq.(6) to the diameter values reproduces a non- 
monotonic behavior of the diameter near Tc (fits 2 and 4 



0.45 




0.6 0.7 0.8 0.9 1.0 1.1 1.2 
T 

FIG. 4: Coexistence curve diameters for bulk LJ fluid with 
L ~ 12 (T (solid circles), bulk LJ fluid with L 20 cr (open 
squares), LJ fluid in pore with H = 12 cr and i « 17 ct (open 
circles) and LJ fluid in pore with H = 12 cr and i ~ 34 a (solid 
squares). The bulk critical densities, reported in Ref. 46] (di- 
amond) and Ref. l28jl (star) are shown for comparison. The 
linear fits 1 and 5 and flts 2 and 6 to Eq.(6) (Table III) are 
shown by dashed lines and solid lines, respectively. 



TABLE IIL Values of the parameters obtained from flts of 
Eq.(6) for the diameter pd to the data for the bulk and con- 
flned LJ fluid, represents the mean-square deviations with 
a confidence limit of 0.95 normalized to a dispersion of the 
simulated data of 0.001. 



N Pc Bl Bi_c 




comments 


bulk LJ fluid, L ^ 


12 a 




1 0.3156(4) 0.2277(15) 


1.23 


highest T = 1.11 


2 0.3269(7) 0.647(32) -0.403(30) 1.15 


aU data 


bulk LJ fluid, L « 


20 a 




3 0.3182(11) 0.2140(67) 


5.2 


aU data 


4 0.3131(35) 0.544(221) -0.30(20) 


4.7 


aU data 


confined LJ fluid, H = 12 


(7, L 


^ 17 a 


5 0.2407(9) 0.2764(41) 


5.6 


highest T=1.09 


6 0.2473(26) 0.492(13) -0.21(12) 


10 


aU data 



in Table III, the former fit is also shown by a solid line in 
Fig. 4). The value pc = 0.3269(7), obtained for the system 
with L « 12 g, noticeably exceeds the critical densities, 
reported in iM El ElEl. The value of the critical 
density, obtained from the GEMC simulations, could be 
underestimated due to the different numbers of particles 
in the two coexisting phases l^^. With increasing 
system size this effect becomes negligible and could shift 
the estimated critical density to higher values. However, 
our simulations show the opposite trend: increasing the 
system size to L « 20 cr causes a decrease of the critical 
density to Pc = 0.3182(11) in the case of a linear fit (fit 3 
in Table III) or even to Pc = 0.3131(35) by fitting Eq.(6) 
(fit 4 in Table III). So, an increase of the system size 
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FIG. 6: Density profiles of the LJ liquid in a slitlike pore with 
H = 12 (7 in equilibrium with a coexisting vapor phase. The 
temperature increases from the top (T — 0.60) to the bottom 
(T = 1.13). 



FIG. 5: Liquid-vapor coexistence curve and diameter of the 
LJ fluid in the slitlike pore H = 12 a, L « 17 a (solid circles) 
and H = 12 (T, L « 34 a (open circles). The bulk coexistence 
curve (line) is reproduced using Eq.(4) with set 4 of the fitting 
parameters from Table II and a linear diameter (fit 1 in Table 
III). 



from i«12cr to 20 cr shifts the value of the criti- 
cal density closer to the values pc = 0.316(3) and pc 
= 0.3197(4) [2^, obtained by finite-size scaling methods. 
Note also, that the mixed-field scaling method, used in 
[28L li^ , overestimates the value of pc of the square- well 
fluid due to the neglect of the Yang- Yang anomaly ^3 • 



2. Pore coexistence curve 

The average density of the LJ fluid confined in a pore 
was calculated assuming that the pore volume is equal to 
L^H. In fact, the fluid- wall interaction is equal to zero at 
the distance 0.86 a from the pore wall and so in an op- 
erational approach this interval could be divided equally 
between the volumes of the fluid and the solid. As a 
result, the real volume accessible to the fluid molecules 
should be reduced by a factor 1.077 and the average pore 
density should increase accordingly. Throughout the pa- 
per we use the average density in the pore without this 
correction on the accessible volume. 

The simulated coexistence densities and the diameter 
for the LJ fluid in the slitlike pore with weakly attractive 
walls of width H = 12 ct are shown in Fig. 5. The data 
points obtained from the simulations using the pore with 
lateral size L « 17 ct and the pore with lateral size L « 34 
a coincide, except at the 3 highest temperatures, where 
the density of the liquid phase becomes lower in the large 
system (see Fig. 5). The critical temperature of the pore 
coexistence curve is depressed due to the effect of the 



conflnement. The vapor phase in the pore has a larger 
density than in the bulk and this effect increases with 
temperature. At low temperatures this may be under- 
stood as a result of adsorption at the weakly attractive 
substrate. At high temperatures this effect is the result 
of a higher coexistence pressure due to the conflnement 
in a pore with weakly attractive substrate 30] . The den- 
sity of the liquid phase in the pore is shifted signiflcantly 
to lower values due to a depletion of the density near the 
weakly attractive surface |3j] and this effects the liquid 
density much stronger than the increasing pressure. 

The pore critical temperature is located between the 
highest temperature, where a two-phase coexistence was 
obtained (T = 1.13) and the lowest temperature, where 
the two phases become identical in the GEMC simula- 
tions (T — 1.15). In a slitlike conflnement the 3D critical 
behavior is distorted by a temperature crossover towards 
two-dimensionality when approaching the pore critical 
temperature and by the spatial heterogeneity of the fluid 
density at all temperatures due to the influence of the 
surface. There are no theoretical equations which give 
a description of the pore coexistence curve. Fits of the 
Eq. (4) to the pore coexistence curve show, that the criti- 
cal temperature varies slightly depending on whether the 
critical exponent /3 is fixed to the 2D or 3D value (see Ta- 
ble II). We determine for the pore critical temperature 
the value T^ = 1.145(2). 

The diameter of the pore coexistence curve is shown in 
Fig. 4. A linear temperature dependence is observed in a 
wide range as in the bulk case (see dashed lines) . The fits 
of the diameter to Eq.(6) are given in Table III (fits 5 and 
6). Deviations of the high-temperature points upwards 
from the linear dependence and convergence to the value 
Pc ~ 0.255 can be noticed. In the larger simulated system 
with L w 34 (T this effect disappears (see Fig. 4). 
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FIG. 7: Density profiles of tlie LJ vapor in tiie slitlike pore 
H = 12 (T in equilibrium with the coexisting liquid phase. 
Temperature increases from the bottom (T — 0.80) to the 
top (T =1.13). 



3. Density profiles 

The density of the coexisting liquid phase in the pore 
is significantly lower than the coexisting hquid density 
in the bulk (Fig. 5). This originates from a depletion of 
the liquid density near the weakly attractive substrate 
and can be studied by an analysis of the density profiles 
where r is the distance from the pore center. At 
low temperatures the density profiles in the liquid phase 
show strong oscillations due to packing effects, which are 
still pronounced even near the pore center (Fig. 6, upper 
curves). Such strong density oscillations prevent observ- 
ing a density depletion at the temperatures T = 0.60 
to T = 0.80. With increasing temperature the density 
profiles become more smooth. However, first and second 
shell oscillations are still to be seen even at the highest 
temperature T — 1.13 (Fig. 6, lowest curve). Visual ex- 
amination of the liquid profiles in Fig. 6 evidences that 
the density near the surface decreases faster with tem- 
perature than the density in the pore interior. 

The temperature evolution of the density profiles of 
the vapor phase is shown in Fig. 7. The vapor density 
maximum near the wall remain pronounced in the whole 
temperature range studied. Its position coincides with 
the first density maximum in the liquid phase profiles and 
corresponds to the minimum of the fluid-wall potential at 
a distance a from the pore wall (see Eq.(3)). In the pore 
interior the vapor density varies smoothly and remains 
flat up to T = 1.07. At higher temperatures the profile of 
the vapor becomes convex upwards and a second density 
maximum can be seen at the highest temperatures due 
to the appearance of packing effects with increasing fluid 
density (Fig. 7). 

To study the surface critical behavior, we define the 
local order parameter Ap(r,r) and the local diameter 
Pd{r,T) of the fluid near a surface similarly to the bulk 
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FIG. 8: Profiles of the order parameter for the LJ fluid in the 
slitlike pore with H = 12 o". Temperature increases from the 
top (T = 0.60) to the bottom (T = 1.13). 
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FIG. 9: Scaling of the order parameter profiles for seven tem- 
peratures from T — 1.04 to 1.10. The dashed line is the func- 
tion Ap{Az)/Aptuik = 1 - exp(-(Az - 0.09 o-)/S,ot''") with u 
= 0.63, ^0 = 0.30 (7. For definition of Az see text. 



case 

Ap{r,T)^ipi{r,T)-p,{r,T))/2 (7) 
Pd{r,T) ^ {pi{r,T) + p,{r,T))/2 (8) 

The proflles of the order parameter Ap{r,T) are shown 
in Fig. 8. Due to a partial compensation of the density 
oscillations in the liquid and vapor phases, the order pa- 
rameter proflles are smoother than the liquid or vapor 
proflles. 

The proflle of the order parameter in a lattice system is 
free from oscillations near the boundary and may be de- 
scribed by smooth fimctions. At distances larger than the 
bulk correlation length ^, the deviation of the order pa- 
rameter from the bulk value decays exponentially, when 
moving away from the surface [51|. Near the surface (at 
distances Az << ^) the order parameter profile should 
obey an algebraic law. Due to the density oscillations of 
the fluid near the structureless surface the algebraic law 
could be observed only in close proximity to the critical 
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point, when the bulk correlation length significantly ex- 
ceeds the range of density oscillations. In a confined fluid 
the approach to the bulk critical temperature is limited 
by the pore critical temperature and therefore the alge- 
braic region may not be detected. That is why we fitted 
the order parameter profiles to a simple exponential law 



Ap(Az, 



1 — exp 



Az-\ 



, (9) 



where Apf,Mife(T) is the order parameter of the bulk fluid; 
^_ is the bulk liquid orrelation length along the coexis- 
tence curve. Taking into account, that distances closer 
then about 0.5 tr to the wall are not accessible to the cen- 
ters of particles, we define Az as the distance to a parallel 
plane 0.5 a inside the pore. As this definition contains 
some ambiguity, we introduced in Eq.(9) the adjustable 
parameter A. 

First, we examine whether A/3(Az,t) really follows the 
universal critical behavior, described by Eq.(9). There- 
fore we plot Ap{Az,T) / Aphuik{T) vs the normahzed dis- 
tance Az/^_ with ^_ — (.ciT^'^, V = 0.63. In the tem- 
perature range 1.04 < T < 1.10 a universal master curve 
perfectly describes the order parameter profiles, exclud- 
ing the first oscillation near the wall (see Fig. 9). The 
dashed line in Fig. 9 corresponds to an amplitude of the 
correlation length = 0.30 a and A = 0.09 a. Note, 
that the center of the first fluid layer is situated near the 
potential minimum at Az — 0.5 cr, while the fluid-wall in- 
teraction is repulsive at Az < 0. So, the order parameter 
becomes equal to zero at some distance between the min- 
imum and the zero value of the fluid-wall potential. At 
higher temperatures (T > 1.10, not shown in Fig. 9) the 
order parameter proflles A/3(Az,r) deviate in the pore in- 
terior from the master curve to lower values. This effect 
is caused by the influence of the opposite wall and will be 
discussed below. At lower temperatures (T < 1.04, not 
shown in Fig. 9) the oscillations of the order parameter 
become more pronounced. This limits the temperature 
range where the validity of the universal behavior could 
be analyzed. 

In a second approach, the correlation length was esti- 
mated from fits of Eq.(9) to each profile Ap(Az,r). The 
obtained values of Apb^ik and ^_ are shown in Fig. 10. 
The parameter A fluctuated around zero and never ex- 
ceeded a few hundredths of cr, when the whole profile 
Ap(Az,t) is used for fitting. If only parts of the pro- 
file Ap(Az,r) are used and the parameter A is fixed to 
0, the obtained values of Apbuifc and ^_ varied no more 
than about 5 to 10 %. For illustration, the values 
obtained from fits to parts of the order parameter pro- 
files extending from the wall to various distances Az are 
shown in Fig. 11 for two temperatures. Comparison of the 
values Aphuik, which were obtained from the fitting of 
Eq.(9), with the values of the true bulk order parameter, 
obtained from the bulk coexistence curve (Fig. 10, upper 
panel), shows, that the two sets coincide only at low tem- 
peratures. With increasing temperature the fitted val- 
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FIG. 10: Temperature dependence of the fitting parameters 
Apbuifc (upper panel) and ^_ (lower panel), obtained from 
fits of Eq.(9) (open circles) or Eq.(lO) (open squares) to 
Ap(Az,r). Tlie solid line in tlie upper panel shows the be- 
havior of the bulk order parameter. In the lower panel the 
dashed line has the slope u = 0.63, the dot-dashed line has 
the mean-field slope i/ — 0.5. The solid line in the lower panel 
is the fitting curve ^_ = O-STt"" (1 - O.S9T° '^^)a. 



ues of Aphuik deviate slightly upward from the true bulk 
data. At the highest 3 to 5 temperatures the fitted values 
of Aphuik rapidly turn down, which correlates with the 
deviations of the Ap(Az,T) profiles at these temperatures 
from the master plot, shown in Fig. 9. At the highest tem- 
peratures the temperature dependence of the correlation 
length ^_ agrees with the expected asymptotic critical 
behavior with an amplitude = 0-43 a (Fig. 10, lower 
panel, dashed line), while at lower temperatures it bends 
downward. Note, that the whole data set of ^„ can be 
fitted by an equation with one correction-to-scaling term 
(see solid line at the lower panel of Fig. 10), using an am- 
plitude ^0 = 0.57 a and a comparatively large amplitude 
(0.89) of the correction-to-scaling term. Thus the values 
^0 = 0.43 (7 and = 0.57 cr of the amplitude of the 
correlation length were obtained when fitting the order 
parameter profiles including the first oscillation. Natu- 
rally, this effectively causes an increase of the obtained 
values ^0- An estimate of from the universal master 
curve (Fig. 9) neglecting the first oscillation of the order 
parameter results in the lower value = 0.30 a. As the 
first density oscillation refiects only some non-universal 
microscopic detail of the interaction potential, the latter 
estimate seems to be more reasonable. For comparison. 
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FIG. 11: Variation of the fitting parameter when Eq.(9) 
is fitted to parts of the order parameter profiles Ap(Az,r), 
which extend from the wall to various distance Az from the 
surface (with fixed A = 0). 



the only estimate of the correlation length for a LJ fluid 
(truncated, however, at another radius) from molecular 
dynamic simulations gives an amplitude of the correla- 
tion length ^0 ~ 0.27 ± 0.2 cr|^], using the critical tem- 
perature of the LJ fluid with full potential [s^- However, 
use of the critical temperature of the LJ fluid with trun- 
cated potential .52.] increases the value almost by a 
factor of two. The available experimental studies of no- 
ble gases give for subcritical temperatures values of £n of 
about 20 to 22% of their Van-der-Waals diameters |53l |. 

The values ^_ , which were obtained from the local or- 
der parameter profiles of a fluid in a slitlike geometry, 
differ from the bulk correlation lengths in semi-infinite 
systems due to two main effects. The shift of the liquid- 
vapor phase transition due to the confinement decreases 
whereas the effect of the opposite wall could in- 
crease the "effective" value of the estimated correlation 
length. Thus, a non- monotonous variation of the "effec- 
tive" value of ^_ with the pore width could be expected. 
To study these effects, simulations of the coexisting den- 
sities in pores of larger width should be performed and 
these studies are in progress now. This should provide 
a new way to study bulk correlation lengths of fluids by 
computer simulations. 

Additionally we have fitted the order parameter profile 
by the mean- field equation 'si'l : 

Ap(Az,T) = Apt,uik{T)tanh (^^^^^ . (10) 

These fits give values of the correlation length ^_ which 
are about 1.6 times lower than the values, obtained from 
the fits using Eq.(9) (see Fig. 10). At high temperatures 
the behavior of the correlation length ^_ is compatible 
both with the Ising model {i/ — 0.63, dashed line in fig. 10) 
and with the mean-field prediction {v — 0.5, dot-dashed 
line in Fig. 10). The behavior of Apbuik obtained from 
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FIG. 12: A double logarithmic plot of the local order param- 
eter, averaged over layers of molecular width (see text for the 
details) vs reduced temperature: first (surface) layer (open 
circles), 2-nd to 5-th layers (squares, triangles, diamond, tri- 
angle down, respectively); 6-th layer (pore interior, solid cir- 
cles) . The slopes of the lines correspond to the exponent (3 = 
0.326 (dashed line), /3i = 1 (dot-dashed line), /3i = 0.82 (solid 
lines) . 



the fits using the mean-field Eq.(lO) agrees well with the 
data for the bulk coexistence curve, excluding the tree 
highest temperatures (see Fig. 10, upper panel). 



4-. Thermal and spatial crossover from hulk-like to 
surface-like critical behavior 

Density profiles of the liquid and vapor phases at the 
pore coexistence curve were used to study the tempera- 
ture dependence of the local order parameter Ap(Az,r). 
A double logarithmic plot of the order parameter Api(r), 
averaged over layers of about molecular width a vs the 
reduced temperature r is shown in Fig. 12. The layer be- 
tween the zero energy point of the fluid- wall interaction 
(0.43 a from the wall, Az = -0.07) and the first minimum 
in the density oscillations (l.Scr from the wall, Az = 1) 
was defines as the first (surface) layer. Each subsequent 
layer spreads over 1 a toward the pore center. The order 
parameter in the surface layer A/9i(r) changes practically 
linearly with r (dot-dashed line in Fig. 12). In the second 
layer the behavior of A/92(t) is close to a power law with 
exponent j3i — 0.82 in a wide range of temperatures (see 
squares and solid line in Fig. 12). In the third and subse- 
quent layers a crossover from bulk-like behavior (dashed 
line. Fig. 12) to surface behavior at smaller r, i.e. closer 
to the critical point is observed. 

The crossover temperature could be extracted from 
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FIG. 13: Lines and symbols are the same as in Fig. 12. The 
data for the 3-rd and subsequent layers are shifted vertically. 
The position of the crossover points from bulk-like to surface- 
like behavior in each layer are denoted by stars. 



this plot as the temperature of the crossing point of the 
two straight lines, which represent bulk behavior with the 
exponent (i = 0.326 (dashed lines in Fig. 13) and surface 
critical behavior with the exponent /? = 0.82 (solid lines). 
These crossover points are shown by stars in Fig. 13 for 
the 2-nd to 6-th molecular layers. 

To study the dependence of the corresponding 
crossover temperature r* on the distance from the surface 
Az in more detail, we performed the same analysis for a 
finer grid of layers of 0.2 cr width. The obtained depen- 
dence is shown in Fig. 14. The values of the correlation 
length ^_ , obtained from fits of Eq. (9) to the order pa- 
rameter profiles Ap(Az,T), are also shown in Fig. 14 (right 
axis, solid circles). Comparison of the two dependencies 
shown in Fig. 14 evidences that the crossover from bulk 
to surface critical behavior is governed by the correlation 
length: the crossover in a definite fluid layer occurs, when 
the correlation length ^_ corresponds roughly to half the 
distance from the layer to the solid surface: Az « 2^_. 

The spatial crossover from bulk to surface critical be- 
havior was also studied by using the effective local critical 
exponent /3'^/''^(Az,AT), which was obtained from fitting 
the temperature dependence of Ap(Az,r) by the equa- 
tion 
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FIG. 14: Layer distance Az vs crossover temperature (left 
axis, open symbols) and temperature dependent correlation 
length 5-(t) found from fits of Eq.9 to the profiles Ap(Az,r) 
(right axis, solid circles). 
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Ap{Az,T)^B^ffTP ^'^'''^') (11) 
The value of l3^-^f{Az,AT) essentially depends on the 
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FIG. 15: Spatial crossover of the local effective exponent 
^''^^(Az, AT), determined from Eq.(ll) applied to the local 
order parameter Ap(Az,r) in different temperature intervals 
AT: 0.72 < T < 1.12 (circles, lowest curve), 1.01 < T < 
1.12 (squares, middle curve) and 1.08 < T < 1.12 (triangles, 
upper curve). The critical exponent of the bulk fluid (3 and 
the surface critical exponent /3i of the Ising model are shown 
by dashed and solid horizontal lines, respectively, r is the 
distance from the pore center, Az the distance from the wall. 



temperature interval AT, used in the fit. It is close 
to the value of some true critical exponent if the fit- 
ting temperature interval AT is outside the crossover 
regions. In Fig. 15 we show the local effective critical ex- 
ponent (3''^f{Az,AT) as a function of the distance from 
the pore center for 3 different temperature intervals AT 
used for the fit. f]''^-^ {Az,AT) from averaging over the 
whole studied temperature range (see circles in Fig. 15) 
shows a pronounced crossover from about /3 = 0.38 in the 
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pore center about 1.1 in the first surface layer. Moreover, 
jS'^f f (^Az, AT) shows oscillations, which correlate with the 
density oscillations, observed in the liquid phase at low 
temperatures. In particular, (3'^-f^{Az,AT) achieves the 
first local maximum of about 1.1 at the minimum of 
the fluid-wall potential. If the fitting temperature in- 
terval is reduced to 1.01 < T < 1.12, the oscillations 
of P^-^-f {Az, AT) become less pronounced and the effec- 
tive exponent varies between 1 at the minimum of the 
fluid-wall potential to about 0.5 in the pore center. In 
the narrow high temperature interval 1.08 < T < 1.12 
the spatial crossover of i3''ff{Az,AT) is rather smooth, 
however, in the first surface layer ^ {Az,AT) is still 
higher (about 0.9) than in subsequent layers (see Fig. 15). 
So, the spatial crossover of P^ffiAzAV evidences that 
the localization of molecules near the wall due to the 
fluid-wall potential distorts the surface critical behavior 
toward higher values of /3i. A similar behavior was ob- 
served for the surface layer of water near a hydrophobic 
wall 0. Away from the surface the density oscillations 
disappear at high temperatures and f3i approaches the 
value f3i — 0.82 of the ordinary surface transition of the 
Ising model. 



5. Density profiles of supercritical fluids 

Density profiles of the LJ fiuid confined in the slitlike 
pore were calculated for various average densities along 
the bulk critical isotherm Tc = 1.1876. At low densities 
the attractive fluid-wall potential caused adsorption of 
molecules near the surface and the profile in the pore inte- 
rior is concave down (Fig. 16, a). With increasing density 
the effect of density depletion near a weakly attractive 
substrate develops and provides density profiles which 
are convex up (Fig. 16, b). The largest density gradient at 
T = Tc roughly corresponds to the bulk critical density 
Pc ~ 0.3 a^. Further increase of the average fiuid den- 
sity in the pore results in density oscillations due to the 
packing effect (Fig. 16, c). In the very dense supercritical 
fluid the density oscillations extend through the whole 
pore (Fig. 16, d). 

The density profiles at the bulk critical temperature 
from Fig. 16 were fitted by Eq.(9), using pbuik and ^+ in- 
stead of Aphuik and respectively. The parameter A 
was found close to zero, when the average pore density 
was below 0.15 or above 0.5 a^. In the intermedi- 
ate density range A passes trough a maximum of A w 
0.22 (7, when the pore average density is equal to 0.27 
a^. Note, that a similarly small value of A was obtained 
from the master curve of the order parameter profiles (see 
Fig. 9). The obtained values of the supercritical correla- 
tion length 1^+ are shown in Fig. 17. The maximum of ^-|- 
is observed, when the average pore density is about 0.29 
a^. When the correction factor, which accounts for the 
accessible volume, is applied to the average pore density 
(see section III. 2), this value becomes practically equal 
to the bulk critical density pc « 0.32 a^. Note also, that 
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FIG. 16: Density profiles of LJ fluid confined in slitlike pore 
with H = 12 (T along the bulk critical isotherm Tc — 1.1876 at 
various average pore densities: p = 0.05 a'^ (a), 0.20 < p < 
0.30 cr^ (b), p = 0.50 cr^ (c) and p = 1.00 (d). The profiles 
b, c and d are shifted vertically. The averaged densities are 
shown by dashed lines. 



the corresponding fitting parameter value pbuik ~ 0.42 
in Eq.(9) noticeably exceeds the critical density of the 
bulk LJ fiuid . 

Also, the temperature evolution of the density profile 
at supercritical temperatures was studied along the pore 
critical isochore pc = 0.255 (Fig. 18). All profiles cross 
in a point, where the density is about the bulk critical 
density. The density gradient from the pore surface to 
the interior decreases with increasing temperature and 
achieves a maximum at about the pore critical temper- 
ature T = 1.145. Accordingly, the correlation length, 
obtained from fits of Eq.(9) to the density profiles, also 
achieves its maximum ^-|_ = 1.87 a at this temperature. 
The parameter A varies in the range 0.21 a < X < 0.27 
(7. So, we may conclude, that the value of about 0.2 cr for 
A determines the effective fluid boundary, when the den- 
sity oscillations and the details of the fluid-wall potential 
are weak or negligible. This is the case for the critical 
isotherm in a wide density interval around the critical 
value and for the critical isochore. Besides, neglecting 
oscillations of the order parameter profile also gives a 
similar value of A at the coexistence curve (Fig. 9). 

The temperature dependence of the average and local 
pore density diameters and of the local density along the 
pore critical isochore is shown in Fig. 19. The average 
density diameter as well as the diameters in the surface 
layer and in the pore interior are continuously crossing 
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FIG. 17: Correlation length obtained from the fits of the 
Eq.(9) to the density profiles along the bulk critical isotherm 
Tc = 1.1876 (Fig.16). 
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FIG. 18: Density profiles of LJ fiuid confined in slitlike pore 
with ¥L — 12 a along the pore critical isochore pc = 0.255 cr^ 
at the temperatures T = 1.145, 1.1876, 1.20, 1.25, 1.30, 1.40 
1.50 and 1.60. 



over to the average and local densities at supercritical 
temperatures. The pronounced cusp in Fig. 19 indicates, 
that the anomaly of the pore diameter (fit 6 in Table 
III) originates from the behavior of the fluid in the pore 
interior. 



IV. DISCUSSION 

Computer simulations of a LJ fluid and water 
evidence a highly universal surface critical behavior of 
fluids near weakly attractive surfaces. The long-range 
fluid-wall interaction suppresses a drying transition and 
therefore the order parameter could be studied in a wide 
temperature range. The temperature dependence of the 
local densities of the coexisting liquid and vapor phases 



near the surface could be described by Eq.(l). Near the 
surface this equation is valid in the whole temperature 
range, where the liquid phase exists: from supercooled 
temperatures to close proximity of the liquid-vapor crit- 
ical point (r > 0.05). This surface critical behavior in- 
trudes into the bulk fluid on distances determined by the 
bulk correlation length This means, that at any dis- 
tance from the surface the local order parameter crosses 
over from bulk-like critical behavior with exponent (3 = 
0.326 to surface critical behavior with exponent /3i = 0.8. 
The temperature of this crossover increases with increas- 
ing distance to the surface. The width of the fluid layer 
which follows the surface critical behavior described by 
Eq.(l) is about 2^_ and diverges when approaching the 
critical temperature as t^'^ . 

The estimated value of the surface critical exponent is 
in good agreement with the value /3i = 0.82 obtained for 
Ising magnets at zero surface field (hi = 0) [s^]. The 
largest deviations from the Ising behavior are observed 
in the first surface layer, where the fluid-wall potential 
causes a strong localization of the molecules both in the 
liquid and vapor phases, which does not disappear even 
at supercritical temperatures (Fig. 18). As a result, the 
surface critical exponent, estimated for the first surface 
layer is about 0.9 for water fl^ and 0.9 - 1.0 for the LJ 
fiuid (Figs. 12 and 15). Obviously, this deviation reflects 
the trend toward mean-field behavior with Pi = 1 due 
to the localization of the molecules near the minimum of 
the fiuid-wall potential. Besides, the effective value of the 
local exponent (3'^f^{Az, AT), determined in a wide tem- 
perature range, varies similarly to the density oscillations 
in the liquid phase due to the packing effect (Fig. 15). 

It may be assumed, that the difference Ap between 
the densities of the coexisting phases near the surface 
should follow the behavior ~ t^^ with f3i — 0.82 also near 
strongly attractive surfaces below the wetting tempera- 
ture. The contribution ~ t^^ being present in both co- 
existing phases below the wetting temperature, as given 
in Eq.(l). Its disappearance in the liquid phase above 
the wetting temperature, when a liquid layer appears be- 
tween the vapor and the solid surface, seems to be un- 
likely. So, the liquid density near the surface (both in 
the liquid phase and in the wetting layer) should include 
the positive contribution ^ t^^ . This means, that Eq.(l) 
could be valid also above the wetting temperature for the 
wetting phase. The absence of a vapor phase near the 
surface does not allow to separate symmetric and asym- 
metric contributions in Eq.(l), and so the extraction of 
the contribution ~ t^'^ is complicated by the presence of 
a linear regular term. In order to distinguish two contri- 
butions with close values of the exponents (0.8 and 1 in 
this case), a close approach to the critical point is needed. 
The available simulation methods may provide the high 
accuracy, necessary to solve this problem [H^] . 

A fiuid near an attractive wall in general should be 
mapped onto the normal transition universality class of 
the Ising model. The normal transition assumes that a 
non-zero surface field hi acts on the spins in the surface 
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FIG. 19: Temperature dependence of the average pore density 
diameter (open triangles) and of the local density diameters 
Pd in the surface layer (open circles) and near the pore center 
(open squares). The corresponding average and local densities 
at supercritical temperatures along the pore critical isochore 
p = 0.255 (T^ are shown by the respective solid symbols. The 
behavior of the diameter for the bulk L,] fluid is shown by the 
dashed line. 

layer. If the surface field is absent (hi = 0), the Ising 
magnets show a surface critical behavior of the univer- 
sality class of ordinary transitions. In the case of fluids, 
some strength of the fluid- wall interaction potential could 
provide equal attraction of the liquid and the coexisting 
vapor phase to the surface in the vicinity of the criti- 
cal point. At this unique fluid- wall interaction neither 
an increase nor a depletion of the density is observed 
at the critical temperature js^. At any other fluid- wall 
interaction a preferential adsorption of one component 
(molecules or voids) occurs. The latter case corresponds 
to hi ^ 0, i.e. to the norma? transition in Ising magnets. 

In the case of the normal transition the magnetiza- 
tion in the surface layer mi in linear-response approxi- 
mation contains regular and singular contributions above 
the critical temperature |^ |^ : 

m+ = mic + AiT + ... + A2-aT^~°'- (12) 

Below the critical temperature these contributions ap- 
pear in the magnetizations of both coexisting phases. 
Due to the symmetry of the bulk Ising magnet the two 
phases contain the equal contributions (with same val- 
ues and signs). However, at subcritical temperatures 
the surface magnetization should also include the term 
~ T^^ , which describes the temperature dependence of 
the magnetization at hi = (this term was overlooked in 
Refs-Ulil). This contribution accounts for the missing- 
neighbor effect and, therefore, has opposite signs in the 
two coexisting phases with magnetizations ni{ and m{^. 
So, the surface magnetization in the two phases along the 
coexisting curve should have the following temperature 



dependence: 

m{ = BiT^^ + mic + AiT + ... + ^a-aT^"", (13) 
= -BiT^^^ + mic + AiT + ... + Aa-aT^"". (14) 

In the presence of a non-zero surface field hi the two 
phases of the Ising magnet are no more symmetrical. 
Moreover, the magnetization near the surface does not 
vanish at the critical point. The order parameter of a 
phase transition is a measure for the dissimilarity of the 
two coexisting phases, which should vanish at the critical 
point. In the case of ordinary and special transitions [sij 
both bulk and surface magnetizations are equal to zero 
at the bulk critical temperature and therefore the mag- 
netization is indeed the proper order parameter for these 
transitions. In the case of an extraordinary transition, 
the surface magnetization vanishes at the surface critical 
point, while at the bulk critical point it remains nonzero. 
Strictly speaking, the bulk magnetization serves as an 
order parameter for the bulk phase transition, while the 
surface magnetization is the proper order parameter for 
the surface phase transition, the critical point of which 
is located at higher temperature than Tc. In the case 
of the normal transition there is a single (bulk) phase 
transition and the order parameter should vanish at the 
bulk critical temperature. This means, that the magne- 
tization could not be used as an order parameter for the 
normal transition, as it was done in the available theory 
of the normal transition |^ . 

We propose to use as the order parameter of the phase 
transition in Ising magnet in all cases: 

Am(z,r) = (m^(z,T) -m"(z,T))/2. (15) 

This order parameter is the deviation of the magnetiza- 
tion in each coexisting phase from the average magneti- 
zation at the same temperature: 

md{z, t) = {m\z, t) -f m^\z, r))/2. (16) 

In particular, the generalised order parameter Am(z,T) 
and the diameter m(i(z,r) have the following temperature 
dependence in the surface layer : 

Ami(r) = Sir'^i, (17) 

TOl,d(T) = TOlC + AlT + ... + Aa-aT^-". (18) 

In Ising magnets the order parameter Am(z) is zero 
everywhere at T = Tc at any surface field, excluding the 
case of the extraordinary transition, where Am(z) van- 
ishes at the surface critical point. Note, that the pro- 
posed generalisation (Eq.(15)) is equivalent to the stan- 
dard definition of the order parameter, i.e. the magneti- 
zation, in the bulk case as well as in the cases of ordinary, 
extraordinary and special transitions, because of the sym- 
metry of the coexisting phases (m^ — - m^^). 

The average magnetization defined by Eq.(16) (i.e. di- 
ameter of the coexistence curve) is equal to zero in bulk 
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Ising magnets and also near the surface with hi = 0. In 
the case of a normal transition nirf(z) reflects the response 
of the system to the surface field and describes the prefer- 
ential adsorption of one kind of spins. The profile md(z) 
remains nonzero at T ^ Tc and continuously crosses over 
to the magnetization profiles in the supercritical region, 
showing critical adsorption. 

We conclude, that in the absence of surface transitions 
the surface critical behavior of Ising magnets belongs to 
the universality class of normal transition, described by 
Eqs.(13,14). The ordinary transition, which occurs at hi 
= can be considered as a particular case of the normal 
transition, when md(z) = 0. The proposed definitions of 



the order parameter and diameter (Eqs.(15,16)) allow to 
map the surface critical behavior of fluids onto the normal 
transition in Ising magnets, preserving isomorphism of 
the critical behaviors of bulk fluids and magnets. 
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